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Abstract 

A hybrid parameterization including contributions of hydrodynamics and of ex- 
pectations from the spectral shape observed in p+p collisions is introduced. This 
parameterization can successfully describe identified hadron spectra over a wide 
range of px in Au+Au reactions at y/s NN = 200 GeV for all centralities. The pa- 
rameters of the hydrodynamic source compare reasonably well to other attempts 
to describe the spectra. The description is obtained using one universal suppression 
factor of the hard scattering component independent of pt and hadron species. For 
the fit results obtained the observed nuclear modification factor for the different 
particles converges to a universal suppression behavior for pr > 6 GeV /c. 
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1 Introduction 



One of the most interesting recent observations in ultrarelativistic heavy-ion 
physics was the suppressed yield of moderately high p T neutral pions in central 
Au+Au collisions at y/s NN = 130 GeV with respect to appropriately scaled 
p+p results [1], in contrast to a strong enhancement observed at lower beam 
energies [2]. This was qualitatively supported by the observation of a sup- 
pression in non-identified charged hadron yields [1,3]. Among the theoretical 
studies on the origin of the observed high pt deficit [4,5,6,7,8,9], most are based 
on the early prediction [10] of the so-called jet quenching, i.e. the energy loss 
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of fast partons induced by a hot and dense quark-gluon-plasma. Others invoke 
initial state parton saturation [8] or final state hadronic interactions [9] . While 
these different scenarios are being discussed, the interpretation is complicated 
by the parallel observation of a comparably large proton and antiproton yield 
at intermediate momenta [11], which are apparently not similarly suppressed. 
These observations have been confirmed by the new measurements at y/s NN = 
200 GeV [12,13], which extend out to considerably higher transverse momenta 
for the pion measurement. 

Simultaneously, there are a number of hints from the RHIC experiments that 
equilibration may be achieved already in an early state and that the behav- 
ior of the system might at least partly be described by hydrodynamics, most 
notably the observation of strong elliptic flow [14,15]. Also the hadron momen- 
tum spectra can be described by hydrodynamic calculations and parameteri- 
zations quite well in the low and intermediate momentum range [16,17,18,19]. 
In particular, the large (anti)proton/pion ratio can easily be explained by such 
calculations. While at very high px the influence of hydrodynamic production 
should become negligible, in the intermediate momentum range there will be a 
smooth transition from hydrodynamic behavior to hard scattering. It is there- 
fore of interest to study the implications of hydrodynamic particle production 
for the interpretation of the pion suppression. 

In the present paper I will attempt to describe hadron spectra with a parame- 
terization combining hydrodynamic components and components similar to 
the original particle production in p+p collisions. This is considered a simple 
approximation to the more complicated real situation in heavy ion collisions, 
where, even if hydrodynamics is relevant for a large fraction of the particle 
spectra, there will always be some non-equilibrium contribution e.g. from the 
surface of the collision zone. In reality, part of this would likely be intermedi- 
ate in shape between true p+p spectra and hydrodynamics, but it is beyond 
the scope of this paper to attempt to describe the full non-equilibrium na- 
ture of these collisions. The calculations should nevertheless provide a better 
understanding of the hadron spectra than pure hydrodynamic calculations. 
Furthermore, the inclusion of an explicit hard component will allow to de- 
scribe the spectra out to high pr- 



2 The model 

2. 1 The hydrodynamical parameterization 

The present paper uses an extension of the hydrodynamical parameterization 
of [16], which includes effects of transverse flow and resonance decays. It was 
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originally based on a computer program by Wiedemann and Heinz [20] , which 
in turn builds upon [21]. It parameterizes the particle source at freeze-out as: 



wiMr Nhydr °- MeY {2Jr+1) 

x/ (^-sinh^(O). (1) 

The integral over £ contains a Woods-Saxon spatial profile with a difluseness 
parameter A and the two Bessel functions which originate from the boost 
of a thermal distribution with a transverse collective velocity. The transverse 
expansion is described by a transverse rapidity, which is parameterized as 
Vt(0 — Vf ' £> i- e - depending linearly on the normalized radial coordinate 
£ = r/R. I fix A = 50 which yields a nearly box-like spatial shape of the 
source, as this has been shown to be a reasonable assumption [22,23,16]. A 
fixed upper limit (£ = 4) has to be chosen for the numerical integration, 
for the distribution used in this work this is effectively equivalent to setting 
the upper limit to infinity. A spin degeneracy factor is included, and there 
is an additional arbitrary normalization factor N^ydro, which is the same for 
all particle species and controls the relative strength of the hydrodynamical 
component. While the spectral shape is a.o. determined by the kinetic tem- 
perature Tkm, the normalization for each particle species is readjusted to the 
chemical temperature T c h em assuming that dN/dy at midrapidity scales with 
the temperature as [21]: 

drith V „?/m 2 2m \ ( m 
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The parameterization attempts to describe the production of pions, protons 
and antiprotons and kaons; it includes contributions from the following reso- 
nances: p, K° s , K*, A, S + A, rj, W, f]'. For the description of baryons a bary- 
onic chemical potential \ib is used, while for strange particles an additional 
strangeness suppression factor \ s is introduced. For each hadron species the 
spectrum of hydrodynamically produced particles fiylroiPT) is given as the 
sum of the direct contribution and those of the appropriate resonances: 
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Fig. 1. Transverse momentum spectrum of neutral pions in p+p collisions at 
■y/s = 200 GeV as measured by the PHENIX experiment [25] . The solid line shows 
a fit of equation 4. 

where W is the invariant mass, Y r the rapidity and Mt t the transverse mass 
of the resonance, and F(W, P r ,p) is the appropriate phase space factor. 



2.2 The p+p parameterization 



As suggested by Hagedorn [24], hadron spectra in high energy p+p collisions 
can be described by a QCD-inspired power law with exponential continuation 
at lower momenta. Such a spectral shape can be written as: 

Pt < Pti 

(4) 




Here n and T are adjustable parameters. T is related to the more frequently 
used parameter p by: po = nT — prx- Fig. 1 shows a fit of this function to the 
neutral pion spectrum measured in p+p collisions by the PHENIX experiment 
[25]. pti = 1 GeV/c was chosen. The function provides a perfect fit over the 
full range of the measured data. As fit parameters C = 280.52 mbGeV~ 2 , 
T = 0.226 GeV and n = 9.89 have been obtained. 

These distributions contain a component due to hard scattering processes, 
which is expected to lead to a power law shape dominating at high px- At 
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low Pt production from string fragmentation, a soft process, will be much 
more abundant and therefore dominate there. For the later description of 
spectra in heavy-ion reactions I will attempt to separate the two components, 
although a clear separation directly from the spectra is not possible. However, 
the exponential shape at low pt is believed to be due to the soft component, 
while at very high pt the shape will turn into a pure power law oc p^ n . The 
separation is of course completely arbitrary when looking at the fit to the 
data in Fig. 1 alone. In fact, there the exponential is not even fitted directly 
to the data and is only determined as the differentiable continuation of the 
shape at higher p T . The value of the separation point p T \ = 1 GeV/c is not 
uniquely determined, but it yields a reasonable inverse slope parameter for 
the exponential. 

The separation is necessary for a comparison to heavy ion collisions, because 
the different physical origin of the two components will lead to different scal- 
ing behavior, as the number of hard collisions increases more rapidly with 
increasing mass of the nuclei than the number of strings, i.e. the number of 
participant pairs. I will assume the soft contribution to be described by the 
exponential: 

sW = c ,(5LzM)" exp (M) (5) 



After subtraction of this part the hard component can be written as: 



Kpt) 
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(6) 



where C s = Ch = C for p+p collisions. 

This parameterization with a subtracted exponential may look artificial. One 
should remember, however, that already the original power-law-like shape sug- 
gested by Hagedorn was an effective parameterization which turns into a power 
law for very high p?. The deviation from the power law at lower momenta, 
also to be seen as a variation of the exponent with momentum, was not de- 
rived from underlying principles, but chosen to describe the experimental data. 
The hard component introduced here is equivalent at high pt and takes into 
account partial contributions of a soft component at intermediate pt- 

The pion spectrum is: 

f<? P \pT) = s{p T ) + h{p T ). (7) 
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Fig. 2. a) Assumed relative contribution of the soft and hard components of the pion 
spectrum in p+p. b) Fraction of the total number of charged hadrons belonging to 
different hadron species. The solid line represents the pion fraction, the dotted line 
the kaons and the dashed line the protons and antiprotons. 



This is essentially only a rewriting of equation 4. The main purpose of this 
separation is the possibility to scale the two components independently when 
increasing the system size by going to heavy-ion collisions. In naive pictures 
of particle production the hard component is expected to increase with the 
number of nucleon-nucleon collisions, as cross sections are small and processes 
are incoherent. The soft component, however, will more likely scale with the 
number of participating nucleons as suggested by the wounded nucleon model. 
The prescription used here leads to a continuous transition between the soft 
and hard component as can be seen in Fig. 2a. The two components are of 
similar magnitude for p? ~ 2 GeV/c. 

To describe the spectra of other hadron species further assumptions are neces- 
sary. The only data from similar beam energies and a reasonably large coverage 
in pt are available from the ISR [26]. One of the major observations in this 
paper was, that the py-dependent ratios of different hadron species were to a 
good approximation independent of the beam energy studied. The kaon, pro- 
ton and antiproton spectra have thus been tuned to reproduce the observed 
ratios. This has been achieved by using for kaons: 

/W(p T ) = 0.44 • tanh(0.46 GeV -1 cpr) • s(p T ) + 0.44 • h(p T ), (8) 



for protons: 

, M , s 1.8-tanh(0.25GeV~ 2 c 2 ^) . . . . 

/pp (jPr) = : ~,\ Ar v -i — " <Pt) + 0-33 • h (p T ), 9 

yy 1.3 + 0.4GeV cp T 



and for antiprotons: 
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_ 1.8(0.3 + 0.4 GeV^cpr) • tanh(0.25 GeV~ 2 c 2 ^) 
[PT) ~ 1.3 + 0.4GeV^ T S[PT) 

+ 0.11-% T ). (10) 



The analytic expressions have been chosen because they can conveniently de- 
scribe the shapes observed experimentally, they have no physical significance. 
The parameters have been tuned to describe the experimentally observed ra- 
tios obtained at the ISR (see Fig. 14 in [26]). They lead to particle ratios as 
shown in Fig. 2b. No error analysis has been done for the parameters used in 
equations 8, 9 and 10. Consequently the related additional uncertainty in the 
final fits has not been taken into account. 



3 Fits to experimental spectra 



Fits of the function 

Ax 

' hydro ~ r J pp 



f (x) (PT) = fSL + fS ) (ii) 



have been performed to spectra of identified hadrons as measured by the 
PHENIX experiment in Au+Au collisions at 200 GeV. Preliminary spectra of 
charged pions, kaons, protons and antiprotons [27] and neutral pions [28] have 
been used simultaneously for a given centrality. The centrality classes used 
are summarized in Table 1 together with estimates of the number of binary 

Table 1 

Average number of NN collisions and participant nucleons for centrality classes as 
used by PHENIX [29]. 



Centrality bin < N coU > < N part > 



0-5% 


1065.4 


± 


105.3 


351.4 ± 2.9 


5-10% 


845.4 


± 


82.1 


299.0 ± 3.8 


10-15% 


672.4 


± 


66.8 


253.9 ± 4.3 


15-20% 


532.7 


± 


52.1 


215.3 ± 5.3 


20-30% 


373.8 


± 


39.6 


166.6 ± 5.4 


30-40% 


219.8 


± 


22.6 


114.2 ± 4.4 


40-50% 


120.3 


± 


13.7 


74.4 ± 3.8 


50-60% 


61.0 


± 


9.9 


45.5 ± 3.3 


60-70% 


28.5 


± 


7.6 


25.7 ± 3.8 


70-80% 


12.4 


± 


4.2 


13.4 ± 3.0 
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collisions < N co u > and the number of participants < N part > [29]. Besides 
the parameters of the hydrodynamical parameterization two normalization 
parameters N s = C s /C and N h = Ch/C of the pp parameterization enter. All 
parameters are summarized in Table 2. Four different sets of fits have been 
performed: 

(1) Pure hydro dynamic fits (i.e. N s = N h = 0) using only data for m T — m < 
1 GeV/c 2 . 

(2) Pure hydrodynamics (N s = N h = 0) using the full data sets. 

(3) Pure fits of the pp parameterization (N hydro = 0) using the full data sets. 

(4) Fits of the full hybrid parameterization using the full data sets. 

The hydrodynamic fits (1) to the low momentum region provide excellent fits 
for all centralities. There is some ambiguity in the choice of T kin and ((3 T ) - 
to some extent a large velocity can be compensated by a small temperature, 
and vice versa. To avoid these ambiguities all later fits have been performed 
with setting T^ n = 120 MeV, which provides good agreement at low tut for all 
centralities. One should note that while especially the flow velocity depends 
on the choice of the freeze-out temperature, the results on the scaling of the 
different particle species as discussed below are not strongly affected. 

Table 2 

Fit parameters of the proposed parameterization. Note, that while the average ve- 



locity is given here, this has been calculated from the transverse rapidity parameter 
rjf which was used in the fits. Fit parameter values for central collisions with their 
fit errors are stated in the right column. 


parameter 


remark 


fit results (central) 


Nhydro 


normalization of hydrodynamical contribution 


7.7 ±3.1 




kinetic temperature 


120 MeV 


(Pt) 


average transverse expansion velocity 


0.529 ±0.015 


Tchem 


chemical temperature 


160.8 ±4.3 MeV 




baryonic chemical potential 


0.17 ±0.03 


\s 


strangeness suppression 


0.96 ±0.15 


N s 


normalization of soft component 


124 ± 100 


N h 


normalization of hard component 


173 ± 20 



Examples of these fits to the spectra from peripheral collisions (70 — 80%) are 
shown in Fig. 3. The pion spectra can be well described already by the pure 
pp parameterization as well as the full fit, the pure hydrodynamics fit, how- 
ever, fails to describe the high tail. While for pions all parameterizations 
provide a good description at lower tjit, for kaons, proton and antiprotons 
discrepancies are seen mostly at low momenta. These discrepancies are largest 
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Fig. 3. Transverse momentum spectrum of pions, kaons, protons and antiprotons 
in peripheral Au+Au collisions (70 — 80%) at y/s = 200 GeV as measured by the 
PHENIX experiment [27,28]. The dotted line shows pure hydrodynamics fits (2), 
the dashed line a fit of the pure pp parameterization (3) and the solid line shows a 
full fit (4) of equation 11. 

for the pp parameterization. As I am most interested in the behavior at high 
momenta, this will not be investigated further. One may take it as a hint that 
low momentum particle ratios extracted from pp measurements at ISR [26] do 
not completely describe even peripheral heavy ion collisions at RHIC, while 
the spectral shape of pions appears to be well described over the range inves- 
tigated here. The reasonable description of all particle species at intermediate 
rriT lends some support to the choice of the pp parameterization. The full fit 
combining hydrodynamics and the pp parameterization provides an excellent 
description of all spectra. Of course, the high momentum behavior of the heav- 
ier particles can not be thoroughly tested from the limited momentum range 
alone. 

Similarly, Fig. 4 shows fits to spectra for central collisions (0 — 5%). Again, 
it can be seen that hydrodynamics alone fails to describe the tail of the pion 
spectra, while the other parameterizations do a reasonably good job. However, 
for the heavier particles the failure of the pure pp parameterization is obvious. 
There are large discrepancies for the kaons over the whole momentum range 
and for protons and antiprotons for large momenta. It is clear that the spectra 
for central collisions can not be described neither by hydrodynamics nor a pure 
pp parameterization alone, however the combination of both provides again 
an excellent description. An equally good description can be obtained with 
this parameterization for all centralities. 
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Fig. 4. Transverse momentum spectrum of pions, kaons, protons and antiprotons in 
central Au+Au collisions (0 - 5%) at y/s = 200 GeV as measured by the PHENIX 
experiment [27,28]. The dotted line shows pure hydrodynamics fits (2), the dashed 
line a fit of the pure pp parameterization (3) and the solid line shows a full fit (4) 
of equation 11. 

4 Discussion 



Although the description of spectra by means of hydrodynamic parameteriza- 
tions is not the ultimate aim of this paper, it is still instructive to study the 
behavior of the fit parameters for the different centralities. Fig. 5 shows param- 
eter values for different fits as a function of centrality. The open symbols show 
results from pure hydrodynamic fits. Those have been performed for lower my 
only with the kinetic temperature being either fixed to T^ in = 120 or 140 MeV 
or varying freely, in addition a fit with unconstrained T^ in has been performed 
to the full available ttit range. The filled circles show the same parameters 
for a full fit of equation 11 including the soft and hard components and as- 
suming Tki n = 120 MeV. It is noteworthy that all these fits show essentially 
very similar results. The chemical temperature (Fig. 5a) is T c h em « 160 MeV 
for all centralities, also the baryonic chemical potential (b) appears to be al- 
most independent of centrality with a value of /^s/T cfeem « 0.15 — 0.2. The 
strangeness suppression (c) is strongest for peripheral collisions with X s « 0.65 
increasing to X s ~ 0.9 for central collisions. The transverse flow velocity also 
increases considerably with centrality from {(3t) ~ 0.25 for peripheral reac- 
tions to (Pt) ~ 0.5 in central reactions. The chemical temperature appears to 
be similar, but slightly smaller than values obtained for central collisions of 
Au+Au at 130 GeV [16,30]. Kinetic temperature and flow velocity are similar 
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Fig. 5. Parameters of the hydrodynamical component of fits to Au+Au collisions 
at ^/s = 200 GeV as a function of the number of participants, a) chemical temper- 
ature, b) baryonic chemical potential normalized to the chemical temperature, c) 
strangeness suppression factor, d) average transverse expansion velocity. The open 
symbols show the results from pure hydrodynamic fits under different conditions, 
the filled circles show results using the full parameterization. 

to values obtained in [16] for 130 GeV or in [18] for 200 GeV, when considering 
that the latter analysis did not take the influence of resonances into account. 
The results of the hydrodynamical parameters obtained here are therefore not 
at all extraordinary or astonishing. (For reference the fit values for central 
collisions using the full fit are also given in Table 2.) 

Already in [16] it was argued that the hydrodynamic contribution to the 
hadron yield even beyond pt = 3GeV/c is not negligible. This question can 
be revisited with the help of the refined fits performed here. Fig. 6 shows the 
relative contributions of the hydrodynamical component to the full fits for 
central collisions. Again, it can be seen that this contribution is significant 
for pions beyond pt = 3GeV/c. Moreover, the importance of this contribu- 
tions reaches out to even higher p? for the heavier particles, as is expected 
for a hydrodynamic source. 2 This hydrodynamic component causes very nat- 
urally a different behavior of the particle ratios compared to p+p collisions. 
The proton/pion and antiproton/pion ratios (as displayed in Fig. 7a) show 
pronounced maxima close to a value of 1 between pt of 2 and 4 GeV/c. Con- 
sequently there is a strong minimum in the pion/hadron ratio (see Fig. 7b). 



2 The strong increase of the relative importance of the hydrodynamic contribution 
in protons and antiprotons at low px is due to the suppression of these particles in 
the non-hydro production (see Fig. 2). 
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Fig. 6. Ratio of the hydrodynamical contribution to the total spectrum in central 
Au+Au collisions at y/s = 200 GeV as a function of pr- The ratios for pions, kaons, 
protons and antiprotons are shown separately. 

The kaon/hadron ratio rises much faster than in p+p. 

This will also influence the comparison of spectra from central nuclear colli- 
sions to those from p+p collisions as is customary in the investigations of high 
Pt hadron suppression. One uses the nuclear modification factor: 

u x , , (l/N%>l)d 2 NX A /dydp T 
AA[PT) (N coll )(l/N^)d*NX/dydp T 1 ) 



This can easily be calculated from the parameterization used here. Fig. 8 shows 
Raa for the different hadron species. Raa for pions decreases from about 0.4 
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Fig. 7. Hadron ratios in central Au+Au collisions at y/s = 200 GeV as a function 
of pt- a) Protons/pions and antiprotons /pions, b) charged pions/charged hadrons, 
charged kaons/charged hadrons and proton+antiprotons/charged hadrons. 



12 



< 1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



012345678 

p T (GeV/c) 

Fig. 8. Nuclear modification factor (equation 12) as a function of px for pions, kaons 
and protons+antiprotons. 

at 2 GeV/c to 0.2 at high p T . For heavier particles Raa starts out at higher 
values and reaches the same asymptotic value more slowly. For protons and 
antiprotons there is a peak structure with a small enhancement (Raa > 1) i n 
the intermediate px range. 

Finally one can investigate the relative strength of the soft and especially 
the hard component, N s and Nh as a function of centrality. Fig. 9a shows 
N s /(0.5 x N part ) as a function of N part . Without taking into account any 
hydrodynamic production this normalized strength of the soft contribution 
increases significantly with increasing centrality, while for the full parameter- 
ization the ratio stays close or below a value of 1. The difference between the 
two estimates may just be seen as the importance of the added hydrodynamic 
contribution, which very naturally appears to be negligible for the most pe- 
ripheral collisions and increases strongly with centrality. A similar difference 
is seen when normalizing to N co u (Fig. 9b), however, now all values are smaller 
than 1. Even more interesting is the evolution of the hard component. Fig. 9d 
shows Nh/N co u as a function of N part . This can effectively be seen as an av- 
erage value of Raa over the pt range of the hard component. Qualitatively 
the trend of both parameterizations is similar, there are, however, significant 
quantitative differences. For the pure pp parameterization the apparent aver- 
age suppression changes from 0.8 for peripheral to 0.4 for central reactions. 
When accounting for the hydrodynamic contribution, the ratio changes more 
drastically from no suppression (= 1) to a value of 0.2. Within this model, 
these are the true suppression values of the hard component independent of 
species and px- 

The normalizations of the soft and hard component are well defined from their 
fit to pp data, the hydrodynamic component has an overall arbitrary normal- 
ization related e.g. to the unkown source volume. Once the normalization is 
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Npart N part 

Fig. 9. Strength of the soft and hard component as a function of the number of 
participants. The open crosses show the pure pp parameterization, the filled circles 
the full fit including the pp parameterization and hydrodynamic contributions, a) 
AV (0.5 * N part ), b) Ns/Ncou, c) AV (0.5 * AT part ), d) N h /N coU 

fixed for a given centrality, it is still instructive to compare it to those for 
other centralities. Fig. 10 shows N hydro /N part normalized to one for the most 
central collisions as a function of the the number of participants. While N part 
varies by a factor of almost 30 over the centrality range investigated, the ratio 
appears to be constant within errors, i.e. the size of the hydrodynamic system 
is to a large extent determined from the collision geometry. This is another 
hint for the consistency of the model used. Within this parameterization one 
can deal simultaneously with a suppression of the hard scattering component 
and a build up of hydrodynamic behavior. It would of course be of great inter- 
est to study explicitly the dynamics of the interactions between jets and the 
hydrodynamic system as e.g. performed in [31] - this is however far beyond 
the scope of this paper. 



5 Summary 

Fits of a hybrid parameterization containing both hydrodynamics and ele- 
mentary soft and hard contributions to hadron spectra measured in Au+Au 
collisions at y/s NN = 200 GeV [27,28] have been performed. The input dis- 
tributions have been tuned to describe the spectra of neutral pions in p+p 
collisions at yfs = 200 GeV [25] and the particle ratios as measured at the 
ISR [26] . Keeping the parameters of the hydrodynamic source and the relative 
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Fig. 10. Strength of the hydro dynamic component per participant as a function of 
the number of participants. The ratio is normalized to one for the most central 
collisions (see Table 2. 

normalizations of the soft and hard components as free parameters an excel- 
lent description of hadron spectra for all centralities can be obtained. The 
parameters obtained for the hydrodynamic source are reasonable: The chem- 
ical temperature of T chem m 160 MeV and the baryonic chemical potential of 
IJB/T c hem ~ 0.15 — 0.2 are almost independent of centrality. For peripheral 
collisions the strangeness suppression factor is A s ~ 0.65 and the flow velocity 
is (Pr) ~ 0.25, while for central collisions I obtain \ s 0.9 and (f3 T ) 0.5. 

The hydrodynamic source contributes a significant fraction of hadrons at inter- 
mediate pt reaching out to at least px = 3 GeV/ c for pions and pr = 5 GeV/ c 
for protons and antiprotons, thereby explaining the large baryon-to-pion ratio 
in central Au+Au collisions. It also results in a very different scaling behavior 
for the different species, even if one assumes one universal suppression factor 
of the hard component for all particle species. The suppression would only be 
visible for proton and antiprotons for pt > 4GeV/c. The universal suppres- 
sion factor extracted is smaller than 0.2 for central collisions, the merging of 
Raa for different species can be seen for px > 6 GeV/c. 
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